SPINVERT: A program for refinement of 
paramagnetic diff'use scattering data 



Joseph A. M. Paddison^^^, J. Ross Stewart^, and Andrew 
L. Goodwin^ * 

^ Department of Chemistry, Inorganic Chemistry Laboratory, University of 
Oxford, South Parks Road, Oxford OXl 3QR, U.K. 

E-maih andrew . goodwinachem . ox . ac . uk 

^ ISIS FaciUty, Rutherford Appleton Laboratory, Harwell Science and Innovation 
Campus, Didcot, Oxfordshire OXll OQX, U.K. 

Abstract. We present a program (SPINVERT; http://spinvert.chem.ox.ac.uk) 
for refinement of magnetic diffuse scattering data for frustrated magnets, spin 
liquids, spin glasses, and other magnetically disordered materials. The approach 
uses reverse Monte Carlo refinement to fit a large configuration of spins to 
experimental powder neutron diffraction data. Despite fitting to spherically- 
averaged data, this approach allows the recovery of the three-dimensional 
magnetic diffuse scattering pattern and the spin-pair correlation function. We 
illustrate the use of the SPINVERT program with two case studies. First we 
use simulated powder data for the canonical Heisenberg kagome model to discuss 
the sensitivity of SPINVERT refinement to both pairwise and higher-order spin 
correlations. The effect of limited experimental data on the results is also 
considered. Second, we re-analyse published experimental data on the frustrated 
system Yo.5Cao.5BaCo407. The results from SPINVERT refinement indicate 
similarities between Yo.5Cao.5BaCo407 and its parent compound YBaCo407, 
which were overlooked in previous analysis using powder data. 
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1. Introduction 

The intense current interest in magnetic materials in which the magnetic moments 
("spins") do not form a periodic arrangement has been fuelled by the observation 
of strongly-correlated states, including the discovery of emergent magnetic charges in 
spin-ice materials [HE] , the observation of excitations with fractional quantum number 
in quantum spin liquids [3], and the presence of chiral order in otherwise-degenerate 
states [1]. States such as these — lacking long-range magnetic order, but exhibiting 
distinctive forms of local order — have been termed "cooperative paramagnets" [S]. 
The "gold standard" experiment for the study of spin correlations is neutron scattering 
on a large single crystal sample. Such measurements provide a direct measurement 
of the Fourier transform of the three-dimensional spin correlation function /(Q). 
In cooperative paramagnets, neutron scattering experiments often reveal highly- 
structured diffuse features in reciprocal space — such as the well-known "pinch points" 
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of spin ice — which can be analysed to characterise the spin correlations in real-space 

mm- 

Unfortunately, however, large single-crystal samples are not available for many 
interesting materials. This is particularly — but far from exclusively — a problem 
for newly-synthesised compounds. Consequently, neutron scattering measurements 
must often be performed on powder (polycrystalline) samples. In a powder neutron 
scattering experiment, the three-dimensional scattering pattern is collapsed onto a 
single axis, Q. For magnetic materials which show conventional long-range spin order, 
this may not actually present a significant problem, since the established technique 
of magnetic Rietveld refinement can (in favourable cases) provide an unambiguous 
solution of the magnetic structure based on only the powder-averaged positions 
and intensities of the magnetic Bragg peaks [7j. In this way, symmetry arguments 
effectively replace the information which is lost as a result of powder averaging. In 
cooperative paramagnets, by contrast, the lack of long-range spin order means that 
there are no magnetic Bragg peaks: the powder average of the scattering pattern shows 
no sharp features, merely a few diffuse "humps" . As a consequence, there has been no 
obvious analogue of magnetic Rietveld refinement which can determine the determine 
the three-dimensional spin correlation function from powder diffuse scattering data. 

In this paper we introduce a new suite of programs, SPINVERT, which allows 
determination of the three-dimensional spin correlation function by fitting powder 
diffuse scattering data. Our program employs reverse Monte Carlo refinement [HI 
to fit a large configuration of spin vectors to experimental powder data. The positions 
of spins are fixed at their crystallographic sites throughout the refinement, while their 
orientations are refined in order to fit the data. The RMC algorithm itself is well- 
established, having previously been used to study states as diverse as spin glasses 
[U, ordered magnetic structures HU], and flux- line lattices in superconductors [TT]. 
The algorithm is entirely analogous to a direct Monte Carlo simulation, apart from 
one important difference: the function which is minimised during the refinement is 
not an energy term defined by a spin Hamiltonian, but rather the sum of squared 
residuals which quantifies the level of disagreement between the fit and experimental 
data. Thus, RMC is not a technique for modelling magnetic interactions based on a 
spin Hamiltonian, but rather a technique for refinement of spin correlations based on 
experimental data [12] . 

It is clear even from this brief overview that the RMC approach is a simple — even 
naive — one. However, for cooperative paramagnets, it is also remarkably powerful. 
In our previous work |13| . we investigated the extent to which three-dimensional 
information could be recovered from powder diffuse scattering data in the following 
way. First, powder diffuse scattering patterns were simulated for a number of models 
of frustrated magnetism. These powder "data" were then fitted using the RMC 
approach. From these refined spin configurations, we calculated the three-dimensional 
scattering pattern /(Q), and compared the patterns obtained from RMC with the 
exact results. We found that, for each model we considered, almost all the features of 
the single-crystal scattering patterns were reproduced by fitting the powder patterns. 
Consequently, the spin configurations obtained by RMC correctly reproduce the 
three-dimensional spin correlations of the starting model: they represent one of the 
degenerate "solutions" of the paramagnetic structure. 

Why is the RMC approach so successful? To answer this question, we compare 
it to alternative "model-independent" approaches. Most such approaches essentially 
involve fitting a simple form for the radial spin correlation function to the powder 
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data. In these methods, the connectivity of the crystal structure is not considered; 
instead, either the magnitude of the spin correlations is fitted for a one-dimensional 
Hst of atomic separations (see, e.g., HHHS]), or the diffuse features are fitted to a 
broad peak-shape function to estimate the spin correlation length (see, e.g., [TS]). 
However, the crystal structure plays no less significant a role in paramagnets than in 
ordered magnets: indeed, in a paramagnet the crystal structure completely determines 
the symmetry and periodicity of the diffuse scattering pattern jl7) . Consequently, 
knowledge of the crystal structure imposes very significant constraints on the form 
of both the single-crystal and powder scattering. Unlike other model-independent 
approaches — but in common with approaches using a spin Hamiltonian — RMC 
refinement uses knowledge of the crystal structure to constrain the spin correlations. 

Our paper is organised as follows. First, we introduce our implementation of RMC 
in detail, and summarise the key equations underpinning the SPINVERT program. 
Next, we discuss possible ways in which the resulting spin configurations may be 
analysed, drawing connections with previous work using RMC and other techniques. 
Finally, we present two case studies intended to illustrate the use of the SPINVERT 
program. In the first case study, we fit simulated data for the Heisenberg model on the 
kagome lattice [THl [T^], and investigate the effects of limited Q-range and statistical 
errors on the results which are obtained. In our second case study, we present a new 
analysis of existing experimental powder data on Yo.5Cao.5BaCo407 (from Ref. [TS]) 
in order to clarify the nature of the paramagnetic spin correlations in this material. 
The results we obtain from SPINVERT analysis indicate some close similarities 
between Yo.5Cao.5BaCo407 and its parent compound YBaCo407, that may have 
been overlooked in previous analyses [TS] . We conclude with a discussion of the general 
advantages and disadvantages of the RMC approach, and some perspectives for future 
work. 

2. Theoretical background 

Reverse Monte Carlo method 

The general RMC method has been described in detail elsewhere [131 HO]; here we 
summarise our specific implementation of magnetic RMC. First, a supercell of the 
crystallographic unit cell is generated. The supercell usually contains several thousand 
atoms, and periodic boundary conditions are used to avoid edge effects. A classical 
spin vector with random orientation is assigned to each atom, and the goodness-of-fit 
to experimental data is calculated: 



Here, I{Q) is the powder-averaged magnetic scattering intensity j21j . superscript 
calc and expt denotes calculated and experimental values, iy{Q) is an experimental 
uncertainty, and W is an empirical weighting factor. A spin is then chosen at random 
from the supercell and rotated by a small amount. This is done by choosing a unit 
vector s with random orientation and forming the new spin vector 



where < A < 1 is the maximum spin move length. In practice, we have found that 
the final results obtained in RMC fitting are almost independent of the value of A; a 




(1) 



I new 



Si + As 
|S, + As|' 



(2) 



SPINVERT paramagnetic structure refinement 



4 



value A — 0.2 is used by default. After a move is proposed, the change in goodness- 
of-fit is calculated, and the proposed move is either accepted or rejected according to 
the Metropolis criteria (35]. This process is repeated until no ftirther reduction in 
is observed. In this way, spins are iteratively rotated to form a configuration for which 
spin correlations are consistent with experimental data. 

Throughout this work, we treat spins as having unit length, with the magnitude 
of the spin (i.e., the length of the effective magnetic moment) absorbed into an overall 
intensity scale factor. When fitting experimental data, it is almost always necessary 
to allow the value of this scale factor s to refine in order to fit the data. The best-fit 
value of s can be calculated after each proposed move by minimising with respect 
to s, which yields [23] 

Eq [/calc(Q)W(Q)]/KQ)]' 
EQ[^calc(Q)]VKQ)]' 

If the data are placed on an absolute scale, the effective magnetic moment is then 
determined by the relation 

y.'' =g^S{S + l) 

(4) 

where g is the g-factor. It is possible also to refine a flat or linear-in-Q background in 
addition to the scale factor; the relevant equations are given in Ref. 



Magnetic neutron scattering intensity 



The magnetic scattering cross-section is variously denoted /(Q), 'S'(Q) or dcr/dfi 
Within the quasistatic approximation and for a single type of spin, it is given as |24) 

,2 1 ■ 
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(5) 



Here, f{Q) is the magnetic form factor, Q is the scattering vector (wavevector transfer) 
and r.i the position of spin . The lack of divergence of the magnetic field means that 
the neutrons only "see" the component of the spin perpendicular to the scattering 
vector, 



Si = S, - [(Q • SOQ] /Q". 
The proportionality constant C is given by 

C= ^^""^ ^ 



(6) 

(7) 
(8) 

where 7n is neutron magnetic moment in nuclear magnetons and r^ is the classical 
electron radius. Multiplying out the squared modulus in Eq. ([5| allows us to write 
the scattering cross-section in terms of correlations between pairs of spins. 



0.07265 barn, 



COS (Q • Ti 



(9) 



where r^j = Vj — r^. Here, the double sum over spin pairs excludes the self-correlation 
terms, for which i = j. These terms give an incoherent contribution to the scattering. 
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which is the additive factor of ^ in ([9j). In an ideal paramagnet, with zero correlation 
between different spins, this incoherent scattering is all that is measured. In a 
real paramagnet, however, some degree of spin correlation is still present, often to 
temperatures far above the temperature of magnetic ordering. 

In order to calculate the powder scattering cross section, we need an expression 
for the spherical average of ([9]). An exact expression was given in the 1960s by Blech 
and Averbach [21] . but appears to have been used only rarely since then. Instead, 
an approximate form — applicable only to an isotropic paramagnet — has often been 
employed. In this particular case we have 

(S,^.S|) = ^(S,.S,), (10) 

so the spin direction is decoupled from Q. The spherical average over is easily 
performed by writing 

Q • r,y = Qr,-j cos 9. (11) 

The result is just a sine Fourier transform of the spin correlations, analogous to the 
Debye formula or a liquid or a glass 

/isotropic(Q) = gC [M/(Q)]^ + ^ E . S/^^^ . (12) 

It is important to recognise that this expression is exact only for an isotropic 
paramagnet in the limit of a large number of spins. It cannot be used to study 



materials with magnetic anisotropy {e.g., Ising spins). Since (12) is not exact for a 
general spin configuration, in SPINVERT we use the exact expression of Ref. |21) . 
We sketch their derivation here. Taking ^ as the starting point, we define a local 
coordinate system for each pair of spins, The local z axis is directed along the 
vector separating the pair of spins; the local x axis is perpendicular to z and lies the 
plane of Sj; and y is the remaining vector in a right-handed set. This gives for the 
coordinate axes: z — Tij/rij, x = [Si — (Si • z) z] / |Si — (S.^ • z) z|, and y = z x x. 
With this local coordinate system, we apply the identity 

• S| = E f^"' - ^) SrS^ (13) 



a, 13 



where, from Fig. [ij the components a,/3 G [a;,?/, z] of Q are given in spherical 
coordinates by 

Q ~ [Qsinf^cos (p, Q sin 6* sin Q cos 0]. (14) 
Substituting (14) into ([l3|, and the result into gives an expression for / (Q) 



/ (Q, 9, if). The integrals over 9 and ip can then be done by hand, which gives the final 
result 



nQ) = cw{Q)f 





cos Qr 



in which 



A, = SfS^, (16) 

B, j ^2S^S] - S^SJ. (17) 




Figure 1. Local spherical coordinate system used for calculation of the powder- 
avcragcd neutron scattering intensity. 



The sum is taken over all pairs of spins which are separated by radial distances 
< rij < T'max! when periodic boundary conditions are imposed, the maximum radial 
distance r^ax is given by half the length of the shortest side of the RMC supercell 
("box"). 

We end this section by anticipating a common criticism of the RMC approach. 
It is sometimes argued that the use of a large supercell means that the number of 
parameters p to be estimated greatly exceeds the number of data points d. The 
number of parameters is often counted as p = 2N for N spins (since each spin has two 
rotational degrees of freedom), which indeed is usually greater than n. It is implicit 
here that the number of degrees of freedom of the RMC fit is calculated using the 
standard expression u = n — p. However, this is only applicable when the function 
which is fitted to the data is linear in the parameters. Since the magnetic scattering 
cross-section depends on the dot product of pairs of spins, this is not the case, and the 
effective number of parameters is not well-defined. Moreover, the quantity of interest 
from Monte Carlo simulations — both direct and reverse — is also not the orientations 
of individual spins, but rather the correlations between spin-pairs. Consequently, the 
effectiveness of RMC is not well assessed by such arguments, and is better judged by 
its performance in real-world applications. 

3. Performing SPINVERT refinements: some guidelines 

In order to perform a refinement using SPINVERT only a few essential pieces of 
information are required. First, we need four pieces of experimental information about 
the system being studied: 

(i) Experimental powder diffuse scattering data. The input data should contain only 
magnetic diffuse scattering. Ideally the data will have be obtained using neutron 
polarisation analysis to isolate the magnetic signal |26j . However, data for which 
the magnetic signal has been obtained by subtraction of a very high-temperature 
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dataset from a lower-temperature dataset can also used. We emphasise that the 
quality of the results ultimately obtained is entirely dictated by the quality of the 
data; some guidelines in this respect are given in Section [5j 

(ii) Crystallographic details: lattice parameters and fractional coordinates of atoms 
in the crystallographic unit cell. The SPINVERT program requires that the 
specified unit cell has orthogonal axes. For hexagonal or rhombohedral systems 
it is necessary to define a larger orthogonal unit cell for use in SPINVERT. It is 
further assumed that all atoms in the unit cell have identical magnetic properties 
(magnetic moments and magnetic form factors), though all atoms need not be 
crystallographically equivalent. The atom positions are not all required to be 
occupied by a spin, so it is possible to refine data for systems with substitutional 
(but not positional) disorder. 

(iii) Magnetic form factor. An accurate parameterisation of the magnetic form factor 
is likely to prove important in SPINVERT refinement, probably to a greater 
extent than in conventional Rietveld refinement |27) . If the magnetic form factor 
is not well approximated by tabulated functions, it may be worth performing an 
accurate measurement to determine it before proceeding. 

(iv) Single-ion anisotropy. Spin anisotropy may be fixed by giving the direction of 
the easy axis for each atom in the unit cell (for Ising spins), or the direction 
perpendicular to the easy plane (for XY spins). Even if Heisenberg spins 
are specified, it is still possible for the refinement to produce anisotropic spin 
configurations, should this be required to fit the data. However, there are almost 
certain to be many more ways of fitting the data in which spins do not point 
along their anisotropy axes, so SPINVERT is likely to underestimate the extent 
of magnetic anisotropy jl3]. Therefore, if other measurements indicate that spins 
behave as essentially pure Ising or XY variables, we suggest constraining the spin 
anisotropy accordingly. 

There are also four main options concerning the SPINVERT refinement itself: 

(i) Weight. The weight, W in ([T]), is equivalent to inverse temperature in a direct 
MC simulation, and determines the proportion of "bad" moves that are accepted 
despite leading to an increase in — that is, how closely the data as a whole are 
fitted. Its value can be freely specified by the user. In practice, the value of W 
should be chosen such that the fit reproduces the real features of the data without 
over-fitting to the experimental errors. A good starting point is often to choose 
W such that between 25% and 75% of the "bad" moves are accepted when the 
refinement has converged. 

(ii) Box size. The size of the spin configuration, given by the number of unit cells 
along each Cartesian direction which make up the RMC supercell. The maximum 
correlation length which can be modelled is given by one half the length the 
shortest side of the box, so a larger box is needed to reproduce sharper features 
in I{Q). However, once a sufficiently large box has been identified, there is 
nothing to be gained by increasing the box size further — doing so merely gives the 
refinement additional degrees of freedom to vary the long-range spin correlations, 
which it may do by fitting to high-frequency components (e.^., noise) in the 
experimental data. 

(iii) Maximum simulation time. In general, refinements should be run until no 
reduction is seen in x^- The number of moves required depends on the data. 
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the system and the refinement parameters. In practice, ^10'' moves per spin is 
often sufficient. 

(iv) Intensity scale and background corrections. When refining experimental data it is 
almost always necessary to treat the overall intensity scale as a fitting parameter. 
This is the case when the data have not been placed on an absolute intensity 
scale, and/or the magnitude of the magnetic moment is not known precisely. 
In some cases, the minimum in can be quite flat as a function of the scale. 
Then the uncertainty in the scale may be the main source of error in quantities 
derived from the refinement (e.g., spin correlation length). The uncertainty in 
derived quantities may be determined, for an assumed uncertainty in the scale, 
by performing separate refinements with the scale factor fixed at (say) ±5% of its 
refined value. It is also possible to refine a flat or linear-in-Q background term; 
however, this should not be done unless necessary, since it introduces an extra 
parameter into the fit. 

The key output from SPINVERT is a spin configuration — a file containing a 
list of spin vectors, with their corresponding positions in the RMC supercell. An 
individual spin configuration represents a snapshot picture of a possible arrangement 
of spins which is consistent with experimental data. In this sense, a spin configuration 
represents a magnetic structure "solution". A direct consequence of the absence 
of long-range magnetic order is that this magnetic structure solution cannot be 
unique: running the reverse Monte Carlo program again will produce a different 
spin configuration, which may be thought of as a second snapshot taken a long time 
after the first. In fact, it is generally important to run RMC simulations several 
times and average calculated quantities, so that a degree of statistical confidence 
can be obtained in the results. In the following section, we consider how the spin 
configurations obtained from RMC can be mined for information about the system 
properties. 

4. Analysis of paramagnetic spin configurations 

Once RMC refinements have been performed, the next — and most interesting — stage 
of the analysis can proceed: unravelling the actual physics which is contained in 
the spin configurations. Since paramagnetic spin configurations are not unique or 
periodic, there is no way to describe them that is at once comprehensive and general 
while using few parameters. This point is equally applicable to disordered spin 
configurations obtained from other techniques (such as direct Monte Carlo simulation) , 
and reflects a general problem of how to characterise disordered states [551 US] • To 
obtain the fullest understanding possible, it is often useful to calculate quantities 
in both real and reciprocal space: the real-space picture identifies the local spin 
correlations, whereas the reciprocal-space picture identifies periodicities and long- 
range components to the correlations [30 . To facilitate this kind of analysis, the 
associated program SPINCORREL plots real-space spin correlation functions, and 
the program SPINDIFF calculates the three-dimensional magnetic scattering pattern 
/(Q). 

Real-space: spin correlation functions 

It is tempting to try examining directly the orientations of spins within the 
configuration, but the lack of periodicity means that meaningful patterns are often 
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hard to discern. Instead, it is better to plot spin-pair correlation functions, which 
provide a measure of how some property of spin-pairs varies as a function of the 
separation of the pair of spins. In general, 

1 ^ 

(/(S(0),S(r))) = -5]/(S„S,), (18) 

i 

where /(Si,Sj) is some function of a central spin and a spin Sj at vector r 
away from it, and an average is performed over all spins in the configuration as 
centres. The powder diffraction pattern is essentially a direct measurement of the 
Fourier transform of the radial spin correlation function (S(0) • S(r)) = (cos(y9), which 
measures the average scalar product of pairs of spins separated by distance r. For 
a spin configuration, this average is performed by taking each spin S.^ in turn, and 
then averaging the spin correlations over its Zij neighbours, Sj, which coordinate it 
at distance r: 

{SiO)-S{r))^—\2 E S'-S„ (19) 
n(r) ^-^ ^-^ 

where n(r) is given by 

N 

n(r)=^Z,,(r), (20) 

with Zij{r) the number of spins j which coordinate a central spin i at distance r. 

There are many other correlation functions which have been invented to 
demonstrate various aspects of structures; some examples are considered in the case 
studies which follow. In structures involving triangular motifs it may be useful to 
plot the correlation function of the vector or scalar chirality, defined in Ref. |19j . 
or the spin nematic correlation function [Section [S] . In polarised neutron scattering 
experiments on single-crystal samples it is possible to separate the different Cartesian 
components of the scattering function with different crystal orientations |2] ; the spin- 
flip and non-spin-flip scattering cross-sections which are measured are also accessible 
from SPINVERT refinements. In Section M we will consider the three-dimensional 
spin correlation function, the real-space analogue of the three-dimensional scattering 
function /(Q). Finally, we note that SPINVERT refinements provide access not only 
to the configurational average of correlation functions, but also to their distributions. 
For example, (cos-;/') — —0.5 may well indicate 120° correlations, but could also result 
from a bimodal distribution of correlations centred at 90° and 180° . 

Reciprocal space: magnetic scattering functions 

The three-dimensional magnetic scattering pattern /(Q) can reveal important 
properties of frustrated magnets that are not always so clear in real-space — for 
example, the power-law correlations in many frustrated magnets appear as "pinch- 
points" in /(Q). In principle, calculating /(Q) from a spin configuration is 
straightforward using ([5| or (|9]). While both these equations are mathematically 
equivalent, it is much faster to compute /(Q) using ^ (which requires 0{N) 
operations) rather than ^ (which requires 0{N'^) operations). However, choosing 
to calculate ([5| leads to some complications in practice. First, it is not possible to 
allow for periodic boundary conditions. This is obviously problematic for refinements 
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which were performed with periodic boundary conditions appUed. Second, when the 
size of the RMC box is large compared to the spin correlation length of the system, the 
calculated scattering pattern appears very noisy. This seems counterintuitive — since 
using a larger supercell will reduce the error in (S(0) • S(r)) — but the problem is that 
as the supercell size is increased, so too does the number of inter-atomic vectors which 
exceed the spin correlation length. The Fourier transform of these long-wavelength 
components appears as high-frequency noise in the calculated scattering pattern. In 
a real experiment, by contrast, this effect is averaged out by the instrumental Q- 
resolution [5T] . 

In fact, both these observations suggest the same thing: that it is necessary 
to restrict the maximum range of the interatomic vectors considered in ([s]). When 
periodic boundary conditions are applied, the maximum interatomic vector should 
not exceed half the box size in each direction. A smaller cutoff may be chosen to 
correspond to the length-scale at which spin correlations become negligible. The 
question of how to impose such a cutoff was originally addressed in 13 1|, where it 
was pointed out that the noise in calculated scattering patterns could be reduced by 
dividing the RMC supercell into a set of smaller regions called "sub-boxes" or "lots" , 
and averaging the scattering intensity over all sub-boxes. In effect, this procedure 
applies a box filter in real space. The calculation proceeds by calculating the three 
equations in sequence: 

Ar(Q) =^SJ,,„exp(iQ.r,), (21) 

a=l 

where is the position of spin r within a crystallographic unit cell, R is a lattice 
vector giving the origin of the sub-box within the RMC supercell, and the sum runs 
over the Na atoms in the unit cell; 

R+L 

Mr(Q)= AR,(Q)exp(iQ.R'), (22) 

R'=R 

where L is the sub-box size; and finally a sum over all sub-boxes: 

7(Q) = ^|Mk(Q)^ (23) 

R 

Here, we have used the lattice periodicity to write Yi —ra + R. 

The method of sub-boxes greatly reduces the level of noise, and produces 
calculated scattering patterns in good agreement with experimental data. However, 
the original method of ^T] also has some disadvantages. In particular, the origin of 
each sub-box was chosen to be a random lattice vector of the supercell with position R, 
so a quantitatively different answer is obtained each time the calculation is performed. 
This problem can be avoided by using every lattice vector R in the supercell as the 
origin of a sub-box. However, a straightforward implementation of such an algorithm is 
very slow. In SPINDIFF, we accelerate the calculation using a fast blurring algorithm 
often employed in gaming and computer graphics. This takes advantage of the fact 
that adjacent lots contain the same set of spins, apart from at the lot edges, so that 
in one dimension we can write 

M,+i(Q) = M,(Q)-f A,+L+i- A,. (24) 



The key point is that applying ( 24 ) along each of the Cartesian directions in sequence 



generates (22) without redundant calculations. 
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While the method of sub-boxes can improve the appearance of the calculated 
/(Q), it remains important to have sufficient statistical averaging. Calculated patterns 
usually appear noisy if they are obtained using only a few thousand spins; to obtain 
good statistics, it is usually necessary to include '^lO^ spins in the calculation. 
The easiest way of achieving this is to average the /(Q) over many individual spin 
configurations. 

As was the case during the SPINVERT refinement, we make some assumptions 
when calculating /(Q): 

• The program calculates only the magnetic scattering; it does not calculate any 
nuclear scattering. 

• The positions of the spins are fixed at their crystallographic positions. 

• There is no magnetoelastic coupling. 

• All the spins have identical magnetic properties (magnetic form factors and 
magnetic moments). 

Provided these assumptions are satisfied, the program SPINDIFF can be used to 
calculate the magnetic diffuse scattering for any spin configuration, including those 
obtained from direct Monte Carlo simulations, spin density-functional theory, and 
other computational approaches. 



5. Kagome Heisenberg magnet 

As a first case study we consider the Heisenberg nearest-neighbour antiferromagnet on 
the kagome lattice. This canonical frustrated model system was studied in detail in the 
early 1990s, and does not show long-range magnetic order at finite temperature. At low 
temperature, spins on each triangle form a coplanar 120° arrangement, and there is a 
macroscopic degeneracy of ways in which these coplanar units can be arranged |32 1 ll8 j . 
Several real materials seem to realise aspects of this behaviour with either classical or 
quantum spins — for example, hydronium jarosite |331I34| and herbertsmithite [51135) — 
although in practice there is usually some subtle feature which breaks the degeneracy 
of ground states. 

In order to demonstrate the capabilities of SPINVERT, we first simulate powder 
diffuse scattering "data" for this model. We consider the classical spin Hamiltonian 
with nearest-neighbour interactions, 

i7 = -^J^S,.S„ (25) 

at a low temperature T/ J = 0.01. Direct Monte Carlo simulations for this Hamiltonian 
were performed for spin configurations of size 3600 spins (comprising ten stacked 
kagome slabs each of size 360 spins) using periodic boundary conditions. The 
simulation time was 10^ proposed MC moves per spin, and 16 independent simulations 
were performed in order to obtain good statistics for the scattering calculations. The 
powder diffuse scattering pattern was then calculated from the MC spin configurations 



using (15), with an arbitrary magnetic form factor (Ho'^"'") applied. In order to 
facilitate comparison between powder and single- crystal diffuse scattering patterns, 
we plot powder scattering patterns as a function of the dimensionless parameter 
Q' — a(5/27r, and label high-symmetry points [hkl) in the single-crystal patterns. 
Powder data were simulated in the range < Q' < 5, and binned at an interval 
AQ' = 0.02. 
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Figure 2. Kagome lattice. The hexagonal unit cell is shown in black, the 
orthorhombic cell used for RMC refinements in green, and local basis vectors 
X used to determine are shown as red arrows. The local x axis for a given atom 
points towards the centre of one of the triangles to which the atom belongs, and 
the local y axis is a perpendicular vector. Note that duo to the point symmetry 
X and —X are equivalent directions, as are y and — y. 

Having obtained these simulated "data" we use them as input data for SPINVERT 
refinements. Since SPINVERT works with orthogonal axes, it is necessary to convert 
the hexagonal unit cell with dimensions to an orthorhombic cell containing twice as 
many atoms; the transformation matrix for this (frequently-encountered) case is given 

by 
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During the following refinements we set the overall intensity scale factor as a parameter 
to be fitted. Although the absolute scale is already known for the simulated data 
considered here, this is very unlikely to be the case with experimental data, so we follow 
the general procedure of allowing a scale factor to refine. We performed SPINVERT 
refinements for 300 moves per spin, after which time no significant improvement was 
observed in the value of x^. The fit-to-data obtained using RMC is shown in Fig.[3]^a), 
showing essentially perfect agreement with the data. In fact, this is usually the case 
for RMC refinements; indeed, if a refinement fails to fit a significant feature in a high- 
quality dataset, this is usually an indication that something is wrong with the starting 
model or the refinement parameters. 

Having obtained spin configurations in agreement with experiment, we consider 
how well these fitted spin configurations reproduce the known features of the kagome 
Heisenberg model. Looking first at reciprocal space, we calculate the single-crystal 
scattering pattern /(Q) from the SPINVERT refinement and compare it with the 
corresponding exact pattern obtained from the direct MC simulations [Fig. |3][b)]. 
Both patterns were calculated using SPINDIFF. Excellent agreement is obtained 
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Figure 3. Left: SPINVERT fit to powder data for the antiferromagnetic 
Heisenberg model on the kagome lattice. Input data are shown as black circles, 
fit red line, and difference (data-fit) as a blue line. Right: Single-crystal 

scattering in the (hkO) reciprocal-space plane. The panel labelled "Exact" is 
calculated from the model used to generate the input data (left image); the panel 
labelled "RMC" is calculated from the SPINVERT refinement of these simulated 
powder data, as discussed in the text. Both panels are on the same intensity scale. 



between the exact result and the SPINVERT prediction, with only a slight blurring 
of the pattern — corresponding to disordering of the spins — notable in the SPINVERT 
pattern. We note that this ability to recover the three-dimensional scattering pattern 
from spherically-averaged data is not limited to this particular model system. In 
our previous work |13| . we considered seven separate systems with different crystal 
structures and magnetic interactions, and found that the three-dimensional scattering 
pattern was recoverable from powder data in each case. 

We now ask how the spin structures obtained by SPINVERT can be characterised 
in real space. In Fig. |4]ja) we plot the radial spin correlation function (S(0) • S(r)) for 
the SPINVERT refinements, and compare it with the exact result. The SPINVERT 
correlation function is nearly identical to the exact result, typically underestimating 
the magnitude of the correlations by a few percent. Thus it appears that powder 
data are very effective at determining the radial spin correlations. While this result is 
perhaps unsurprising, it is not necessarily the case, since I{Q) and (S(0) • S(r)) only 
contain equivalent information in the limit of large Qmax- 

Following on from this, a natural question is whether the RMC refinements might 
afford some sensitivity to higher-order spin correlations. Such correlations are indeed 
present in the kagome Heisenberg AFM model. At low temperatures {T/J < 0.02), 
short-range spin nematic order develops, so that adjacent triangles of spins tend to 
lie in a common plane; the mechanism may involve "order by disorder" [181 119) . 
The "standard" spin correlation function is insensitive to nematic order, so, following 
Ref. [18) . we define a nematic correlation function which is equal to 1 for a coplanar 
state: 

g,(r) = ^((n(0).n(r)f)-i, (26) 

where 

n= ^(Si X S2 + S2 X S3 + S3 X Si) (27) 
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Figure 4. (a) Radial spin correlation function (S(0) ■ S(r)). Exact values (from 
direct Monte Carlo simulations) are shown as black squares, and values obtained 
from the SPINVERT fit to powder data are shown as red circles. The SPINVERT 
values are near-identical to the exact values, (b) Nematic spin correlation function 
gii{r) defined in the text. Exact values are shown as black squares, and values 
obtained from the SPINVERT fit are shown as red circles. The SPINVERT fit 
correctly determines the trend in nematic correlations, but underestimates the 
magnitude by ~50%. 



is the vector normal to the plane formed by the three spins of a triangle. The 
comparison of the SPINVERT gn{r) with the exact result is shown in Fig.Qb). It is 
apparent that SPINVERT consistently underestimates the value of gnir) by around 
50%. Thus it seems that SPINVERT is somewhat less sensitive to higher-order spin 
correlations than to spin-pair correlations. Conversely, since the experimental data 
are only directly sensitive to spin-pair correlations, it is perhaps surprising that RMC 
should show any sensitivity to higher-order correlations at all. 

Finally, we consider the directional dependence of the real-space spin correlations. 
The correlation functions we have plotted so far are radial averages, which take no 
account of the geometry of the kagome lattice. However, this must obscure a great deal 
of information — after all, it is the geometry of the lattice which leads to frustration. 
In order to represent the directional dependence of the correlations, we consider the 
two-dimensional spin correlation function (S(0) • S(r)). This function is seldom shown 
in the literature (although its Fourier transform /(Q) is frequently given), but it can 
reveal surprisingly simple patterns in the paramagnetic structure |36j . For a Bravais 
lattice, (S(0) • S(r)) is straightforward to calculate: one simply calculates the dot 
product of a central spin with the spin at displacement r away, and averages over 
each spin as a centre. For a non-Bravais lattice, such as kagome, the calculation is 
slightly more involved, because different basis atoms can be related to each other 
by rotational as well as translational symmetries. In such cases, it is necessary to 
rotate our coordinate axes when moving between basis atoms such that each atom 
is in an equivalent environment. For the kagome lattice, with a three-atom basis, a 
set of suitable local axes are shown in Fig [2j The directional dependence of the spin 
correlations is shown in Fig. ([5]). The "Exact" spin correlations for the starting model 
[Fig. (l5|, left panel] show a non-trivial dependence on both distance and angle, yet the 
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Figure 5. Two-dimensional spin correlation function (8(0) ■ S(r)) for the 
Heisenberg kagome antiferromagnet. The left panel shows exact values, and the 
right panel shows values obtained from SPINVERT refinement of powder data. 
The X and y axes are the local axes given in Fig. [2] 



result from SPINVERT refinement reproduces all the main features of this dependence 
[Fig. ([5]), left panel]. 

An obvious question which remains is: how well do these results transfer to 
actual experimental data? Compared with the simulated data used above, real data 
are likely to have a somewhat more limited Q-range, and to contain a certain amount 
of statistical noise. We will now consider the effects of these limitations, and make 
some suggestions as to the kind of data which are suitable for SPINVERT refinement. 

We investigate the effect of limited Q-range on the predicted single-crystal 
scattering patterns by identifying two extreme scenarios. First, we attempt a 
refinement omitting all the data except the main diffuse peak, by setting QJ^ax = 2.5 
[Fig. |6][a)]. Then we consider the inverse situation, omitting all the data below 
Qmin = 2.5 but keeping the subsequent peaks [Fig. [6]jb)]. As anticipated, in both 
cases the reconstruction of the single-crystal scattering is of lower quality than for the 
full refinement. What is surprising is that both refinements nevertheless reproduce 
the main features of the diffuse scattering. Moreover, the overall quality of both 
reconstructions is similar, even though the input data do not overlap in Q. This 
behaviour can be understood by recognising that the diffuse scattering pattern is 
periodic in reciprocal space, and this periodicity is dictated by the crystal structure 
|17| . Consequently, the powder-averaged scattering from subsequent Brillouin zones 
contains information about the first Brillouin zone, and vice versa. In effect, every 
measured value of Q represents an average over a different selection of points in the 
first Brillouin zone, so that increasing Qmax places increasingly stronger restrictions 
on /(Q). 

Does this mean it is best to measure data to the highest possible Qmax? In 
practice, the answer is probably only a qualified "yes". An exception is that the 
behaviour of the single-crystal scattering in the limit Q — >■ is uniquely determined by 
the powder scattering as Q — 0. This limit is generally less experimentally accessible 
for larger values of Qmax- As a result, if there is insufficient information in the data as a 
whole to constrain the low-Q behaviour, there is a tendency for SPINVERT refinement 
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Figure 6. (a) Fit to powder data and predicted single-crystal scattering pattern 
for input powder data with limited Qj^^x = 2.5. Input data are shown as black 
circles, fit as a a red line, and difference (data-fit) as a blue line, (b) Fit to powder 
data and predicted single-crystal scattering pattern for input powder data with 
Q'miu ~ 2-5. (c) Fit to powder data and predicted single-crystal scattering pattern 
for input powder data containing statistical noise, as described in the text. 



to refine inaccurate values for I{Q — 0). An example of this behaviour can be seen 
in Fig.[6]3. Another consideration is that, as the scattering becomes weaker at higher 
Q, the experimental errors become proportionately larger, and fitting to errors in the 
high-Q data may produce noise even at low-Q in the predicted single-crystal patterns. 
We may therefore ask how low the statistical errors need to be. To test this, we added 
a small amount of statistical noise to our simulated data and repeated the SPINVERT 
refinement. The noise was simulated by choosing random numbers from a Gaussian 
distribution with standard deviation a — 0.01 (corresponding to an error of 1% in the 
most intense data point). With these (relatively small) errors, we find that the single 
crystal pattern predicted by fitting the powder data is noticeably reduced in quality 
[Fig.|6];c)], although the main features are still recovered. 

We therefore draw two main conclusions regarding the data quality which is 
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desirable for SPINVERT refinement. First, the differences between Fig. [3] and 
Fig. [6]Jc) suggest that it is important to measure beyond just the first peak in the 
scattering pattern. Second, it is essential that the statistical errors in the data are low; 
specifically, they should be low enough that the modulation in the powder scattering 
pattern is well resolved even at high values of Q. While these criteria may be difficult 
to achieve with samples that are small or very highly absorbing (particularly if the 
material in question also has a small S) , in most other cases they should be attainable 
using current-generation diffuse-scattering diffractometers such as D7 (ILL) [57] or 
DNS (Jiilich) [35]. Refinement of experimental data using SPINVERT is therefore 
likely to be both practical and informative — a point we hope to demonstrate in the 
following section. 

6. Legacy data: Yo.5Cao.5BaCo407 

As our final case study we demonstrate one of the key uses we envisage for the 
SPINVERT program: analysis of legacy data. We consider spin correlations in the 
compound Yo.5Cao.5BaCo407 [15]. Based on powder diffuse scattering data, this 
compound was originally suggested to realise a kagome-lattice model similar to that 
studied in Sectionjs] 15j. However, subsequent single-crystal experiments on the parent 
compound YBaCo407 suggested a different picture [31]. Here we attempt to reconcile 
these two sets of results, by using the original powder data of Ref. [H] to predict the 
single-crystal data obtained in Ref. 39J. 

The crystal structure of Yo.5Cao.5BaCo407 (hexagonal PG^mc [^ or trigonal 
P31c [41,) has a unit cell containing eight magnetic Co ions, spread over two 
inequivalent sites (Col and Co2). The Col sites form triangular layers and the Co2 
sites kagome layers, with the two kinds of layers alternating along the c axis, identifying 
the potential for spin frustration [Fig.[7]. The Col atoms are in the Co'^"'" valence state, 
and the Co2 atoms are mixed- valence Co'^^ and Co^+ ; both Co sites are tetrahedrally 
coordinated by oxygen. The spin correlations in Yo.5Cao.5BaCo407 were first studied 
in detail in Ref. |15| . In that study, polarised-neutron powder diffraction data showed 
the absence of periodic magnetic order to T = 1.2 K. It was considered that this lack 
of magnetic order suggested that the kagome planes were essentially decoupled. This 
was explained by postulating that the Col (triangular lattice) ions are low-spin S = 0, 
rather than high-spin S = 2 as might be expected. The experimental basis for this 
model was that the magnetic moment obtained from fitting the neutron scattering 
data was significantly lower than the calculated value assuming that both Co sites 
were high-spin. However, several subsequent studies have suggested a need for this 
conclusion to be revised. First, an inelastic neutron scattering study [41] shows that 
a high-frequency component to the spin fluctuations is present, which lies beyond the 
dynamic range of the measurements perfumed in Ref. [15| . When this component to 
the scattering is taken into account, the expected moment for a model with both Co 
sites high-spin is recovered. Second, X-ray absorption spectra measurements for the 
parent compound, YBaCo407, indicate that both Co sites are high-spin [32]. It was 
also noted that in symmetry an = state for Co requires two spins to align 
parallel in a ^2 orbital, in disagreement with Hund's first rule. To allow a S = 
state a very large distortion from Td symmetry would be required, whereas the crystal 
structure, as determined by neutron diffraction, shows only a small distortion to C31, 

Perhaps most significant is a study reporting single-crystal neutron scattering 
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Figure 7. Arrangement of magnetic Co ions in Yo.5Cao-5Ba C04O7. The Co2 
(Co2+/Co3+) sites form ka gome layers, sliown in blue. The Col (Co^+) sites 
form triangular layers, shown in red. The structure can be described as a network 
of corner-sharing trigonal bipyramids. 

measurements of the parent compomid, YBaCo407 39,. Unlike Yo.5Cao.5BaCo407, 
oxygen-stoichiometric YBaCo407 undergoes a transition to an antiferromagnetic state 
at Tn — 110 K [43j. Just above Tn, at T — 130 K, broad diffuse features were 
observed in the (hkO) reciprocal-space plane, indicating short-range magnetic order 
within the kagome planes. Most interestingly, much sharper magnetic peaks were 
also observed in the perpendicular scattering plane, identifying the presence of longer- 
range correlations along the c axis. These features were successfully reproduced using 
a spin Hamiltonian in which both Co sites are magnetic: it was found that Col spins 
align ferromagnetically with a long correlation length, whereas the Co2 kagome spins 
remain in a spin-liquid-like state. 

An outstanding question is then: do the spin correlations in Yo.5Cao.5BaCo407 
resemble the model of independent kagome planes of Co2 spins proposed in Ref. |15] , or 
do they follow the model of ferromagnetic correlations between Col spins identified in 
Ref. [53] for YBaCo407? The lack of single-crystal samples of Yo.5Cao.5BaCo407 
renders a direct answer inaccessible. However, we will show that a reasonably 
conclusive answer can nevertheless be obtained by employing SPINVERT refinement 
of the original powder data of Ref. [TS]. We note the existence of a previous RMC 
study Yo.5Cao.5BaCo407 using different data and refinement software [41]; our results 
are in broad agreement with this study, but go beyond it by addressing the three- 
dimensional nature of the spin correlations. We performed SPINVERT refinements 
using a 10 X 6 X 6 supercell of the orthorhombic unit cell containing 5760 spins 
(crystallographic details are taken from Ref. [10] )• In our refinements, we assume 
both Co sites are magnetic. As a starting point, we also take the magnetic moments 
on both Co sites to be of equal magnitude; refinements were also performed assuming 
relative spin lengths S = 2 (for Co^"*") and S = 3/2 (for Co^+), which gave similar 
results. Refinements were performed for 300 proposed moves per spin, after which time 
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Figure 8. Powder magnetic diffuse scattering data for Yo.5Cao.5Ba C04O7 from 
Scfiweika et al. The temperature of data collection for the powder data T = 1.2 
K. Experimental data are shown as black circles, SPINVERT fit as a red line, and 
difference (data-fit) as blue circles. 

no significant improvement was observed in and 16 independent refinements were 
performed in order to obtain good statistics for scattering calculations. An intensity 
scale factor was allowed to refine to fit the data, from which we estimate 5' ~ 1.1 
(taking g = 2), in good agreement with the results of [U [T3]. The fit-to-data we 
obtain is shown in Fig. [Sj All the features of the data are captured very well, apart 
from a slight peak broadening in the fit due to the finite box size. 

Having obtained spin configurations in agreement with the experimental data, we 
proceed to characterise the nature of the spin correlations. First we use the SPINDIFF 
program to calculate the single-crystal scattering pattern from the SPINVERT spin 
configurations. This allows a direct comparison between our prediction of the single- 
crystal scattering for Yo.5Cao.5BaCo407 and the experimental single- crystal data 
for YBaCo407 from Ref. [39] . This comparison is shown in Fig. [9j Strikingly, the 
experimental data and the SPINVERT prediction are almost identical. Particularly 
significant is that the sharp reflections at (^02) are present in both the SPINVERT 
refinement and the single-crystal data. We note that the SPINVERT refinements are 
performed starting from entirely random spin configurations, so these sharp features 
cannot appear by chance: they must be required by the input powder data. Therefore, 
a model of decoupled kagome planes with conventional antiferromagnetic correlations 
(as considered in Section [5]) is ruled out. Indeed, there are qualitative differences 
between the powder diffuse scattering for such a model [Fig. [s] and the experimental 
powder data shown in Fig. |8] 

Following Ref. [39) we plot the spin correlations along two directions in real 
space. In Fig. [loj we show correlations along hexagonal [100]-type directions within 
the kagome plane (blue circles), and along the [001] directions (red squares). The 
correlations along [100] are antiferromagnetic and relatively weak. The correlations 
along [001] are always ferromagnetic, and are significantly stronger, with exponential 
correlation length ^ « 13 A. In these essential features, the correlation functions we 




Figure 9. Comparison of predicted magnetic diffuse scattering for 
Yo.5Cao.5BaCo407 with experimental data for YBaCo407. Two reciprocal-space 
planes are shown, (hkO) [left-hand image] and (Okl) [right-hand image]. For each 
plane, the left panel shows the predicted single-crystal scattering based on fitting 
the T = 1.2 K powder data of Schweika et al. for Yo.5Cao.5BaCo407 1151 . and 
the right panel shows the experimental single-crystal data of Manuel et al. for 
YBaCo407 at T = 130K [39). 
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Figure 10. Spin correlation function for Yo.5Cao.5BaCo407 calculated from 
SPINVERT. Red squares show ferromagnetic correlations along the [001] real- 
space direction (c-axis). Blue circles show antiferromagnetic correlations along 
(lOO)-type directions within the kagome planes. Filled points are values obtained 
from refinements assuming magnetic moments of equal length for both Co sites; 
hollow points assume 5 = 2 for Co^+ and 5 = 3/2 for Co^+. 



have shown agree closely with the model for YBaCo407 (Fig. 3 of Ref. [22)- We note 
that a quantitative comparison is not possible with the correlation functions shown 
in Ref. [32]: their results were obtained from direct Monte Carlo simulations at an 
arbitrary low temperature, and hence will show much stronger short-range order than 
we observe for Yg 5Cao.5BaCo407. However, all the qualitative features, both in real 
and reciprocal space, are in excellent agreement. 

In conclusion, we have found that the spin correlations in Yo.5Cao.5BaCo407 
closely resemble those in YBaCo407. The degree of resemblance is remarkable, 
given that there is a difference in the temperature of data collection of two orders 
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of magnitude between the single-crystal data of Ref. [39] and the powder data 
of Ref. [TS]. This result may be explained by the substitution of Ca for Y in 
Yo.5Cao.5BaCo407, which introduces disorder into the exchange pathways, in effect 
reducing the strength of the magnetic interactions. This disorder may also help 
explain why Yo.5Cao.5BaCo407 fails to order magnetically, although the mechanism 
by which order is prevented is not yet understood. Perhaps most importantly, 
the observation that the key features of single-crystal data can be reproduced by 
fitting powder data for a real system presents several opportunities. In the family 
of compounds based on YBaCo407, it was recently suggested that chiral order may 
be present in the paramagnetic phase |44j : the appropriate correlation functions are 
readily calculable using RMC techniques, so this intriguing possibility could be in 
principle be investigated using SPINVERT. More generally, we would argue that the 
initial uncertainty about the nature of the spin correlations in Yo.5Cao.5BaCo407 
may have arisen partially from a desire to impose a model onto experimental data, 
rather than vice versa. The essential advantage of the reverse Monte Carlo approach 
is that the refinement is entirely data-driven, so there is no requirement to think in 
terms of paradigms of frustrated magnetism {e.g., independent kagome sheets) when 
performing initial data analysis. Reality may be more complex than the models we 
imagine; the RMC approach recognises this, and presents us with results which are 
essentially model- independent. 

Conclusions 

In this paper we have detailed the SPINVERT program for refinement of paramagnetic 
powder scattering data using a reverse Monte Carlo approach. This method is 
capable of reconstructing the three-dimensional magnetic neutron scattering pattern 
from spherically-averaged data [13]. In addition, all the spin correlation functions 
of the system are accessed. The RMC technique is primarily sensitive to spin- 
pair correlations; however, in favourable cases we recover information about higher- 
order spin correlations. We have not considered here whether information on single- 
ion anisotropy can be obtained, preferring to treat the anisotropy as a structural 
constraint. However, we anticipate that the technique may also have a limited degree 
of sensitivity to single-ion effects in systems with relatively low crystal symmetry. 

A key point about RMC is that it is entirely independent of a spin Hamiltonian. 
This may be an advantage or a disadvantage: an advantage, because it is not necessary 
to assume a form of the Hamiltonian to model the spin correlations; a disadvantage, 
because it does not produce such a microscopic model as output. Compared to other 
model-independent techniques for analysis of diffuse scattering (such as simple curve 
fitting), we believe that the RMC approach is superior in both the quantity and 
accuracy of information it provides. Most importantly, the crystal structure is a key 
ingredient in RMC refinements. Knowledge of the crystal structure places strong 
constraints on the form of the magnetic correlations, and it is this extra information 
allows the three-dimensional scattering function /(Q) and spin correlation function 
(S(0) • S(r)) to be reconstructed. In fact, the key point here may not relate directly 
to RMC at all: it is simply that the magnetic powder diffraction data — together 
with knowledge of the the crystal structure — contain a great deal of information. It 
seems very likely that other techniques could be developed to access this information 
effectively; for example, it is possible to fit interaction parameters directly to powder 
data within a mean- field theory model (see, e.g., [45tl46]). 
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There are also certain disadvantages to the RMC approach. Most importantly, 
RMC tends to produces the most disordered spin configuration compatible with 
experimental data [12] . In practice, we have not found the extent to which disorder 
is overestimated to be too great, provided the spin anisotropy is specified for Ising 
or XY systems [IS]- However, it is not possible to know in advance the extent to 
which disorder will be created for a given dataset. It is important always to remember 
that the success of RMC depends entirely on the quality of the data with which 
it is provided. The data should show the whole of the first broad peak as well as 
some of the subsidiary peaks, and should ideally have excellent statistics. The lowest 
possible background is desirable, and for this reason polarised-neutron diffraction data 
have significant advantages [37] • We also anticipate that RMC may be less successful 
for crystal structures which are complex — in the sense that they contain many radial 
separations between pairs of atoms which nearly overlap — since the effect of the limited 
real-space resolution imposed by Qmax will then become more significant. 

We hope that the SPINVERT program will be useful in several situations. One 
possible use is the re-analysis of legacy data, as exemplified in Section [6] More 
importantly, when large single crystal samples are currently unavailable, single crystal- 
like information can be obtained using SPINVERT by fitting high-quality powder data. 
A few topical examples should illustrate the potential for discovery. For example, the 
S* = 1/2 kagome system kapellasite was recently synthesised [17], with measurements 
of the bulk properties and inelastic neutron scattering suggesting a quantum spin- 
liquid phase [48]. In the absence of single-crystal samples, SPINVERT refinement of 
powder data would allow detailed characterisation of these spin-liquid correlations. 
Second, "the metallic spin ice" system Pr2lr2 07 [31] exhibits a range of seemingly 
unique properties, including a spontaneous Hall effect in the paramagnetic phase 
suggestive of chiral order [T. Only small single crystals (~lmm'^) are reported 
[50] , so SPINVERT refinement of powder diffraction data would allow a valuable 
comparison of the predicted /(Q) with microscopic models [51]. Furthermore, the 
ability of SPINVERT to use a structural model containing non-magnetic impurities 
would allow studies of canonical frustrated spin glasses such as SrCr9a;Gai2-92:Oi9 
(SCGO) (see, e.g., [S1[S3]) and Ba2Sn2ZnGaio-7^Cr7^022 (BSZCGO) [53]: a model 
of the structural disorder could be produced using, e.g., pair distribution function 
techniques [55], and this structural model used as a starting-point for SPINVERT 
refinement of magnetic diffraction data. In this way, the impact of the structural 
disorder on the magnetic correlations could be compared with theoretical models 
[5^ . Finally, in some cases it may even be advantageous to perform experiments 
on powder samples: for example, for high-throughput studies, for highly-absorbing 
samples [57], or when only small single crystals are available. Consequently, we believe 
that the SPINVERT code — available at http://spinvert.chem.ox.ac.uk — will enable 
useful information to be obtained for a much wider range of cooperative paramagnets 
than is presently the case. 
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